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ABSTRACT 

We investigate the observational signatures of a straightforward evolutionary scenario for protoplan- 
etary disks in which the disk mass of small (<J 50/im) particles decreases homologously with time, but 
the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral 
indicators of the existence of disks, their structure, and mass evolution that may be tested with the 
upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide 
range of masses, 10~ 8 M Q < Mdisk < 10 _1 Mq. Our Monte Carlo radiative equilibrium techniques enable 
us to explore the wide range of optical depths of these disks and incorporate multiple, anisotropic dust 
scattering. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already 
known from millimeter and radio observations. As the disk mass decreases, the excess emission of the 
disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. 
At near-infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of 
disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses 
(K — L) provide a robust means of detecting disks in star clusters down to A/disk ~ 1O~ 7 M0, while the 
far-IR excess probes the disk mass, the caveat being that large inner disk holes can decrease the near-IR 
disk emission. 

Varying other disk parameters (outer radius, flaring, dust size distribution) alter the SED quantita- 
tively, but do not change our general conclusions on the evolution of SEDs and colors with the mass 
of small particles in the disk. Reducing the disk mass results in a clear progression in color-color 
diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus- 
Auriga sources in the context of decreasing disk mass. Different viewing angles yield degeneracies in 
the color-mass relationship, but highly inclined disks are very faint and red and are readily identified in 
color-magnitude diagrams. 

Subject headings: radiative transfer — scattering — accretion, accretion disks — ISM: dust, extinction 
— stars: pre-main-sequence 



1. INTRODUCTION 

Protoplanetary disks are primarily detected through 
their signature infrared excess emission relative to a stel- 
lar photosphere (e.g., Mendoza 1968; Rydgren, Strom, & 
Strom 1976; Cohen & Kuhi 1979; Rucinski 1985). The IR 
excesses arise from the reprocessing of stellar radiation and 
the liberation of accretion luminosity within the disk (e.g., 
Lynden-Bell & Pringle 1974; Adams & Shu 1986; Adams, 
Lada, & Shu 1987; Kenyon & Hartmann 1987; Adams, 
Emerson, & Fuller 1990; Lada & Adams 1992; Hillenbrand 
et al. 1992). At other wavelengths, disk masses and ve- 
locity structures have been measured with radio and mm 
interferometers (Beckwith et al. 1990; Beckwith & Sar- 
gent 1993; Koerner, Sargent, & Beckwith 1993; Koerner 
& Sargent 1995; Dutery et al. 1996; Wilner & Lay 2000) 
and more recently, high resolution HST observations and 
ground based adaptive optics and speckle imaging tech- 
niques have imaged disks via scattered light (O'Dell, Wen, 
& Hu 1993; Burrows et al. 1996; Stapelfeldt et al. 1998; 
Padgett et al. 1999; Krist et al. 2000; Cotera et al. 2001; 
Grady et al. 2000; Roddier et al. 1996; Koresko 1998). 



The disks that are detected range in mass, structure, 
and formation mechanisms. Massive, optically thick, flar- 
ing circumstellar disks, A/disk ^ 10~ 3 M Q , are a natural 
product of the star formation process (e.g., Shu, Adams, 
& Lizano 1987) and are detected around pre-main se- 
quence Classical T-Tauri stars. Low mass "debris disks" 
(Afdisk < 10~ 5 M Q ) are detected around older main se- 
quence stars (Backman & Paresce 1993). The dust that 
gives rise to the IR excesses and scattered light images of 
debris disks is believed to form from collisions of planetes- 
simals within the disks. 

The advent of large IR detector arrays now allows for 
simultaneous multiwavelength observations of large num- 
bers of disks in star clusters (e.g., Haisch, Lada, & Lada 
2000). The upcoming SIR TF mission will greatly increase 
the wavelength coverage for studying clusters, thereby en- 
abling detailed studies of disks over a wide range in age 
and mass. We therefore need to know what spectral indi- 
cators are best suited for detecting disks and determining 
their mass evolution. 

In this paper our working hypothesis is that through 
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time the disk mass decreases but the disk structure and 
stellar luminosity do not change. In our models the dom- 
inant opacity source for A <; lOO^m is small (<; 50^m) 
grains, so another way of stating our model is that we as- 
sume the small particle mass traces the total mass of gas + 
dust. The small grain mass may decrease by the accretion 
of disk material or coagulation and growth of large grains 
and rocks while the gas is either depleted or accreted onto 
the star keeping the gas/dust ratio constant. The recent 
detection of if 2 in (3 Pictoris' disk and the derived gas/dust 
<~ 100 is consistent with our assumption (Thi et al. 2001). 

We investigate the observational signatures of this sim- 
ple model for disk evolution and present radiative equilib- 
rium models for a range of disk masses and viewing angles. 
Model spectra are compared with semi-analytic approxi- 
mations for optically thick flat and flared disks and pre- 
dictions are made as to what spectral regions are most 
sensitive to the evolution of disk mass and at what wave- 
lengths disk detections are robust. We perform simula- 
tions for a fiducial disk structure and dust size distribution 
that fits the SED of the edge-on disk system HH30 IRS. 
(Wood et al. 2001). § 2 describes the ingredients of our 
models, § 3 presents our SED simulations and shows how 
SIRTF colors are sensitive to disk mass. § 4 presents mod- 
els that explore variations in parameters other than disk 
mass and discusses alternative disk evolution models. § 5 
compares our models with color-color diagrams of Taurus- 
Auriga sources and we summarize our findings in § 6. 

2. MODELS 
2.1. Radiative Equilibrium Calculation 

Theoretical temperature distributions and emergent 
spectra have been calculated for T Tauri disks using a 
variety of techniques. These include approximations for 
opaque flat disks (Adams & Shu 1986; Adams, Lada, & 
Shu 1987; Lada & Adams 1992) and extensions of this 
approach to flared disks (Kenyon & Hartmann 1987; Chi- 
ang & Goldreich 1997, 1999); two dimensional radiation 
transfer techniques (Efstathiou & Rowan-Robinson 1991); 
diffusion approximations (e.g., Sonnhalter, Preibisch, & 
Yorkc 1995; Boss & Yorke 1996); disks which are modeled 
as "spherical sectors" (Men'shchikov & Henning 1997); 
vertical structure calculations performed by dividing the 
disk into plane parallel annuli (Calvet et al. 1992; Bell 
et al. 1997; Bell 1999; D'Alessio et al. 1998, 1999); 
and Monte Carlo techniques (Wolf, Henning, & Stecklum 
1999). Common approximations in many of the "tradi- 
tional" codes are to conduct the radiation transfer in only 
one direction in the disk and assume the dust grains scatter 
radiation isotropically. Also many techniques are limited 
to the study of optically thick and hence massive disks. 

Monte Carlo techniques are straightforward to adapt to 
any geometry, mass, and can accurately include polar- 
ization and multiple, anisotropic scattering. We use the 
Monte Carlo radiative equilibrium technique of Bjorkman 
& Wood (2001) which conserves energy exactly. For very 
optically thick disks we use our Monte Carlo technique for 
the upper layers of the disk (i.e., the disk "atmosphere") 
and a diffusion approximation for the densest midplane 
regions. For the most massive disks we simulate, typically 
less than 1% of the flux requires a diffusion treatment. 
More details of our adaptation of the Bjorkman & Wood 



(2001) technique for simulating T Tauri disks will be pre- 
sented in Bjorkman, Wood, & Whitney (2001, in prepa- 
ration). An advantage of using Monte Carlo techniques 
for studying a range of disk masses, is that we are not 
restricted to one-dimensional radiation transfer and can 
therefore simulate SEDs of low mass optically thin disks in 
which radial transport of photons is important. The out- 
put of our code is the disk temperature structure (due to 
heating by stellar photons and accretion luminosity) and 
the emergent SED and polarization spectrum at a range 
of viewing angles. A calculation of the hydrostatic disk 
structure (e.g., Chiang & Goldreich 1997; D'Alessio et al. 
1999) can be included in the Monte Carlo technique, but 
this would require an iterative scheme. At present we have 
not implemented such a scheme and instead perform the 
radiative equilibrium calculation for a fixed disk geometry. 

2.2. Disk Structure 

The determination of disk structure from fitting SEDs 
and scattered light images does not yield a single structure 
and dust size distribution that applies to all disks. Some 
systems are fit with passive flat disks (e.g., Adams et al. 
1987; Adams ct al. 1990; Miyake & Nakagwa 1995), while 
others require flared disks heated by starlight and accre- 
tion luminosity (e.g., Kenyon & Hartmann 1987; Burrows 
et al. 1996; Stapelfeldt et al. 1998). The scattered light 
image of the edge-on disk of HH30 IRS (Burrows et al. 
1996) has for the first time allowed the vertical structure 
of a protoplanetary disk to be studied directly. We there- 
fore adopt the HH30 IRS disk as our fiducial model for our 
SED models. A fixed disk density structure which fits the 
scattered light images (Burrows et al. 1996; Cotera et al. 
2001) and SED (Stapelfeldt & Monetti 1999; Wood et al. 
2001) of HH30 IRS is 

f R \ a 1 
P = Po[-^-J ^-^[z/Kw)] 2 , (1) 

where w is the radial coordinate in the disk midplane and 
the scale height increases with radius, h = h (m/R*) 13 . 
For the HH30 IRS disk we adopt (3 = 1.25, a = 2.25, and 
h = 0M7R+, giving /i(100AU) = 17AU. 

In our simulations the inner edge of the disk is trun- 
cated at the dust destruction radius, J?dust- Assuming 
T* = 4000K, i?* = 2Rq, and circumstellar dust subli- 
mates at 1600K, then i?dust ~ 8i?+. This is larger than 
the dust destruction radius for optically thin dust because 
the reprocessed emission from the disk provides additional 
heating over and above the direct stellar radiation, increas- 
ing the size of the dust destruction zone (a detailed dis- 
cussion of the shape of the dust destruction region will 
be presented in Bjorkman et al. 2001, in preparation). 
In currently popular magnetic accretion models the disk 
is truncated at a radius Rq, which may not be equal to 
-Rdust- If Ro < -Rdust there will be a gas disk extending 
from Rq to i?dust which may give rise to additional IR 
emission. Our models therefore assume that any material 
within i?dust is optically thin, which is a good approxima- 
tion for low mass disks. For high mass disks the gas may 
be optically thick producing larger near-IR excesses. 

2.3. Adopted Circumstellar Dust Properties 
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Recent modeling of HST images (Cotera et al. 2001) 
and the SED of HH30 IRS (Wood et al. 2001) indi- 
cates that the circumstellar dust size distribution extends 
to larger grain radii than typical ISM grains. This is in 
agreement with many other observations indicating grain 
growth within protoplanetary disks (e.g., Bcckwith et al. 
1990; Beckwith & Sargent 1991). This paper primarily in- 
vestigates the effects of disk mass on the SED and adopts 
circumstellar dust properties that reproduce the HH30 IRS 
SED. The dust model (chemical composition, mathemat- 
ical form for the size distribution, calculation of opacity 
and scattering parameters, etc) is described in Wood et 
al. (2001) and we only summarize the main features here. 
Specifically, we adopt, a size distribution 



n(a) da — Cia p cxp (— [a/a c ] q ) da , 



(2) 



with p = 3.5, q = 0.6, a c = 50/im, a m ; n = 0.01/im, and 
flmax = 1mm. The exponential scalelength, o c , yields dust 
particle sizes extending up to and in excess of 50/xm. Fig- 
ure 1 shows the wavelength dependence of the opacity, 
scattering albedo, and Heyney-Greenstein phase function 
asymmetry parameter (Heyney & Greenstein 1941) for this 
size distribution. 

Recent Infrared Space Observatory spectra of Hcrbig 
Ae/Be stars (e.g., Meuss et al. 2000; Chiang et al. 2001; 
van den Ancker et al. 2000) are now allowing the circum- 
stellar dust chemistry to be studied. Because we adopt the 
dust properties of Fig. 1 for our simulations, we have not 
investigated the effects of different chemical compositions 
on the the resulting disk SEDs. How the circumstellar 
chemistry effects the SED is an interesting problem, but it 
is beyond the scope of this paper and we present models 
for different disk masses and the dust properties shown in 
Figure 1. In addition, our models do not include additional 
heating from transiently heated small grains (see the new 
radiation transfer code of Misselt et al. 2001) which is not 
important in the cooler Classical T Tauri stars considered 
in this paper. 

2.4. Energy Sources 

The energy input to the disk is from stellar photons and 
accretion luminosity liberated in the disk. As discussed 
in the previous section, we fix the disk structure for our 
radiation transfer simulations. Given the disk structure 
(Eq. 1), a disk theory determines the accretion rate for a 
given disk mass. Our parameterization of the disk den- 
sity and accretion follows that presented in the review by 
Bjorkman (1997), apart from the term ±y / R /zu. The 

accretion rate and viscosity parameter, M and ayisk, are 
related to the disk parameters by 



M = s/m^ a disk V c Pa h 3 /R* 



(3) 



where the critical velocity V c — ^/GM*/i?+. The flux 
due to viscous disk accretion, GM+M /2R+, is gener- 
ated throughout the disk midplane region according to 
(Shakura & Sunyaev 1973, Lynden-Bcll & Pringle 1974) 



dE _ 3GM+M 
dA dt 47rro 3 



1 - 



(4) 



For low mass disks, the heating due to accretion luminos- 
ity is negligible and stellar irradiation dominates the disk 
heating. In general we choose ayisk = 0.01 (Hartmann et 
al. 1998), but for the most massive disk we simulate this 
results in a very large accretion luminosity, La > 0.8L*. 
For this case ayisk is adjusted so that La < 0.2L* in 
line with recent observational determinations of accretion 
luminosities in classical T Tauri stars (Hartmann et al. 
1998). 

3. MODEL RESULTS: SEDS, COLORS 

The following models use = 4000K, i?„ = 2R Q , 
-Rdisk = 100AU, and a distance of 500pc to the system. 

3.1. SED Evolution with Disk Mass 

Figure 2 shows the effect on the SED of changing the 
disk mass, but keeping the disk structure fixed to that of 
our fiducial model. The massive optically thick disks pro- 
duce SEDs that resemble that of the Kenyon & Hartmann 
(1987) flared disk model, aside from differences (addition 
of scattered light to pole-on views and silicate features) 
due to our inclusion of a finite albedo and non-gray opac- 
ity. The most massive disk has a very large near-IR excess 
due to the large accretion luminosity present in this model. 
The model SEDs display the characteristic features present 
in other simulations: large infrared excess emission, flat 
spectrum sources at intermediate inclinations, and double 
peaked spectra (optical and far-IR peaks) for very large 
inclinations. At large inclinations, the optical peak is due 
to scattered starlight as the dense disk totally obscures the 
star and disk emission at short wavelengths. 

Reducing the disk mass has the most dramatic effect 
at long wavelengths with the SED rapidly declining with 
decreasing mass. At short wavelengths the disk remains 
optically thick to stellar photons over a wide range of mass 
so the near-IR excess is not as sensitive to mass. Figure 4 
presents our models again, with the three panels show- 
ing SEDs for the range of disk masses at a given inclina- 
tion. These results show similar features to the models of 
Men'shchikov & Henning (1997, their Fig. 12) who pre- 
sented SEDs for a range of optical depths in a spherical 
geometry with evacuated bipolar cones. 

The SIRTF sensitivity limits show that at a distance 
of 500pc very low mass disks (M<jisk <; 10~ 6 M Q ) are de- 
tectable out to a wavelength of 70^m. At 25^m SIRTF 
will be sensitive to photospheric flux levels at 500pc, thus 
allowing for the detection of even lower mass disks. Our 
overall conclusion from these simulations is that near-IR 
excesses detect disks while far-IR excesses can be used to 
study their mass. 

3.2. Color Evolution with Disk Mass 

Figure 4 shows the variation of colors with disk mass 
for i < 60°. At long wavelengths we use the simu- 
lated flux at 70pm and 160/im in forming the color and 
have not adopted any particular color system. The col- 
ors are fairly insensitive to inclination for i < 60°, but 
see §3.3 for color-color and color-magnitude diagrams that 
include highly inclined disks. For a passive disk with 
M disk = 1O _1 M , A(K - L) w 0.7, decreasing to a fairly 
constant A(K-L) w 0.4 for 1O" 7 M < M disk < 1O" 3 M . 
Larger K — L colors arise in massive disks where accretion 
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luminosity is included, A(K— L) w 1 for Mdisk = 10 M Q . 
In our models, less massive disks do not sustain large ac- 
cretion rates (see §2.3) and the K — L excess is due to 
reprocessing of starlight. For these disks, if — L is fairly in- 
sensitive to mass as the disk remains optically thick in the 
near-IR for masses Md; s k ^> lO~ 7 Af . Therefore ground 
based near-IR observations are capable of detecting very 
low mass disks. Again, we emphasize that this is the mass 
of small (<; 50/xm) particles which dominate the near-IR 
opacity. 

At longer wavelengths the SED decreases with decreas- 
ing disk mass and this is reflected in the other color indexes 
remaining relatively flat with disk mass until the disk be- 
comes optically thin at the waveband under study. At the 
longest MIPS wavelength, 160/im, the K — 160 color shows 
a clear progression from the most massive to least massive 
disk. 

3.3. Color-Color and Color- Magnitude Diagrams 

Figure 5 shows color-color diagrams for our model disks 
assuming intrinsic stellar colors from Kcnyon & Hartmann 
(1995, Table A5). When we include all disk inclinations 
there are degeneracies in the color-mass parameter space, 
with large viewing angles generally leading to redder col- 
ors. This result differs from other investigations (e.g., 
Lada & Adams 1992; Kenyon, Yi, & Hartmann 1996; 
Meyer, Calvet, & Hillenbrand 1997) which assumed the 
disk emission was proportional to cosi, giving the bluest 
colors for edge-on viewing. More detailed radiation trans- 
fer modeling shows that a simple cos i scaling docs not give 
the correct inclination dependence for flared disk models. 

The effect of changing the disk mass is not very strong 
in the (K — L)/(K — N) diagram for M disk £ 1CT 6 M , 
reflecting that at short wavelengths the disks remain op- 
tically thick over a wide range of mass. The disk masses 
separate out more clearly in the (K — N)/(K — 70) dia- 
gram. Therefore color-color diagrams of large numbers of 
sources that compare near and far infrared colors provide 
a means of determining the range of disk masses within a 
cluster. 

Figure 5 showed that inclination effects lead to a large 
spread in the location of different disk masses in the color- 
color diagrams. This leads to degeneracies in the lo- 
cation of different disk masses. However, Fig. 6 shows 
how color-magnitude diagrams can help in breaking the 
mass-inclination degeneracy. Placing all our disk models 
on color-magnitude diagrams shows that disks viewed at 
large inclinations occupy the lower righthand corners, i.e., 
highly inclined disks are faint and red. 

4. SED AND COLOR DEPENDENCE ON OTHER DISK 
PROPERTIES 

4.1. Flat Disks 

The previous sections presented SEDs for disks which 
have the flaring parameters that fit HH30 IRS images and 
SED. However, other investigations indicate a variety of 
disk structures in Classical T Tauri stars (e.g., Kenyon & 
Hartmann 1987; Miyake & Nakagawa 1995; Chiang et al. 
2001). To investigate the effects of flatter disks on SEDs 
and colors we repeated the simulations of Fig. 2 with a re- 
duced scaleheight h = 0.003.R* giving /i(100AU) = 3AU. 
The smaller scaleheight results in the disk intercepting less 



stellar radiation (e.g., Kenyon & Hartmann 1987) and con- 
sequently smaller excesses and bluer colors than our mod- 
els in § 3. For a given disk mass, the mm flux is unal- 
tered, but the mid to far-IR emission is sensitive to the 
disk structure (Kenyon & Hartmann 1987; Chiang et al. 
2001). The qualitative variation of the SEDs and colors 
with disk mass remains the same. 

4.2. Grain Growth 

In our models the dust size distribution does not change 
as the disk mass decreases. However, small dust grains 
may coagulate to form larger grains and rocks thereby al- 
tering the grain size distribution in the disk (e.g., Beck- 
with, Henning, & Nakagawa 2000). This will lead to 
changes in the wavelength dependence of the opacity. 
D'Alessio et al. (2001) investigated grain growth by keep- 
ing the disk mass constant and increasing the maximum 
grain size (effectively reducing the population of small 
grains) and changing the slope of the power law size dis- 
tribution. For a slope p = 3.5, their SED models show the 
largest effects occur for A ;> 20/im, while forp = 2.5 (which 
puts more mass into larger particles) the SEDs are similar 
to ours in which the total disk mass decreases. SEDs aris- 
ing from either a low mass disk or a massive disk with large 
grains may be distinguished with detailed SED modeling. 

4.3. Disk Radius 

Decreasing the disk radius, but keeping the mass con- 
stant results in a disk that is denser and optically thicker 
than our models of § 3. A smaller disk effectively removes 
cool material at large radii that provides the bulk of the 
long wavelength emission, yielding smaller far-IR and mm 
fluxes (e.g., Beckwith et al. 1990). The other effect of 
squeezing the disk into a smaller volume is an increase in 
the height above the midplane at which stellar photons 
are absorbed. This gives a geometrically thicker dust disk 
that can intercept more stellar photons, raising the near 
and mid-IR emission. The overall effects of a very small 
disk are to increase the wavelength at which the disk be- 
comes optically thin yielding larger K — L colors and a 
slower decline in the K — 180 color with decreasing mass. 
The larger optical depth for small disks will lower the min- 
imum disk mass that may be detected from K — L colors. 
Large radius disks give the opposite of these effects: less 
near and mid-IR, and more far-IR and mm emission. 

4.4. Inner Disk Holes and Non-homologous Disk 
Evolution 

In order to investigate general trends, we have presented 
a simple model for disk evolution of protoplanetary disks 
in which the mass decreases homologously with time. Real 
disks will be more complicated than this, with the disk ra- 
dius growing, changes in the flaring parameters due to dust 
settling, and the possibility of accretion being terminated 
through the opening up of large inner disk holes and gaps, 
resulting in "inside-out" evolution. Current data does not 
provide much support for the existence of large holes in 
protoplanetary disks: if they do not show near-IR excesses, 
mid-IR and mm excesses are usually absent also (Stas- 
sun et al. 2001, Haisch, Lada & Lada 2001a), but there 
are some exceptions notably CM Aur (Kocrner, Sargent, 
& Beckwith 1993). SIRTF should be able to distinguish 
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between low mass protoplanetary disks, whose mass has 
evolved homologously (without creating large inner holes), 
and disks which have evolved by creating large inner holes 
and thus directly test the inside-out disk clearing scenario. 

Debris disks often exhibit large inner holes or a ring-like 
structure (e.g., Koerner et al. 1998; Jayawardhana et al. 
1998; Schneider et al. 1999). However, debris disks are not 
likely to be protoplanetary disks that have evolved large in- 
ner holes. Debris disks are sufficiently old that the dust in 
them is not remnant protoplanetary material but instead is 
created and continually replenished by collisions between 
planetessimals formed previously in the disk. Moreover, 
protoplanetary disks which simply evolve from the inside 
out by clearing large inner holes (but otherwise maintain- 
ing their initial structure) , will still contain large reservoirs 
of material in their outer regions and thus should be con- 
siderably more massive than debris disks. 

5. COMPARISON TO OBSERVATIONS: DISK MASSES IN 
TAURUS-AURIGA 

Observations at mm and radio wavelengths, where the 
disk is optically thin, provide the best probe of disk mass 
(e.g., Beckwith et al. 1990). When such data is not avail- 
able our models show that far-IR, and in certain circum- 
stances even mid-IR data, can provide another means of 
probing disk mass (if the disks are not too massive and op- 
tically thick). For example, Kenyon & Hartmann (1995) 
presented a compilation of optical through far-IR obser- 
vations of Taurus- Auriga sources and noted a pronounced 
gap in the K — N distribution between the bluest Class II 
and the reddest Class III sources. Figure 7 shows the 
(K — L)/(K — N) color-color diagram for the Taurus- 
Auriga sources along with our face-on model colors for 
various disk masses. In the context of our models, the 
gap in the K — N distribution corresponds to disks with 
Miisk <, 1O _6 M0 and indicates that there are few circum- 
stellar disks with masses between 10 -6 and 10 -8 Mq in 
the Taurus-Auriga cloud (of course, this disk mass does 
not include rocks and planets). This, in turn, could be 
interpreted to indicate that once the disk mass falls be- 
low ~ 10 _6 M Q the timescale for clearing the remaining 
material in the inner (<J 5 — 10AU) disk is very rapid. 

6. SUMMARY 

We have investigated the observational signatures of an 
evolutionary model in which the disk mass decreases ho- 
mologously, but the disk structure and stellar parameters 
remain constant with time. Our main conclusions are: 

• Near-IR observations detect disks. Disks remain op- 
tically thick in the near-IR over a wide range of disk 
mass resulting in measurable K — L excesses for disks 
down to Mdi s k ~ 10~ 7 M Q . This corresponds to a 
dust mass of M<j us t ~ 1O -9 M0. 

• Mid-IR observations probe disk structure. Observa- 
tions in the 20yum <J A <; 100/mi range are sensitive 



to the disk flaring parameters and are crucial for de- 
termining the degree of dust settling. 

• Far-IR most sensitive to disk mass. The far-IR emis- 
sion decreases rapidly with disk mass resulting in a 
strong correlation between K — 160 and disk mass. 
Our simulations indicate that at a distance of 500pc, 
SIRTF will be able to detect M disk £ 1O~ 6 M with 
a 500s exposure at 70/im and even lower masses at 
25/im. 

• General trends are reproduced for a range of disk pa- 
rameters. Despite degeneracies in fitting SEDs of 
particular sources, we find that varying the disk pa- 
rameters within generally accepted limits introduces 
a spread in the color-mass relationship, but does not 
affect the general trends. 

The degeneracies inherent in fitting disk models to wave- 
length restricted datasets illustrates that it is necessary to 
have full SED coverage to constrain the disk structure and 
dust size distribution (e.g., Men'shchikov & Henning 1997; 
Chiang et al. 2001; D'Alessio et al. 2001). SIRTF, SMA, 
and ALMA will provide detailed SEDs for nearby sources 
enabling detailed modeling of the disk mass and struc- 
ture. For studies of stellar clusters, such data will not be 
available and mm observations, which are most sensitive 
to disk mass, cannot as yet achieve the required sensitiv- 
ity and resolution to study low mass disks and/or distant 
clusters. Therefore, we must appeal to statistics and de- 
termine average properties of sources in different clusters 
for studying disk evolution. SIRTF will provide broad 
band colors for hundreds of sources within a given cluster, 
from which we can construct median colors and the spread 
around the median for each cluster (e.g., Haisch, Lada, & 
Lada 2001b). From our models it is apparent that NIR ex- 
cesses detect disks, and therefore the absence of a K — L 
excess within a cluster of a given age can determine the 
maximum disk lifetime. By comparing the median colors 
from different clusters with a spread in age (as a whole 
and as functions of position within the cluster), we can 
search for trends as a function of cluster age. Therefore, 
the far-IR probed by SIRTF may allow us to determine 
a timescale for disk evolution — i.e., in our homologous 
evolutionary model we would be able to convert the color- 
mass correlation to one which tracks disk mass in clusters 
as a function of age, analogous to the recent ISO analysis 
of debris disks by Spangler et al. (2001). 
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Fig. 1. — Dust parameters for a grain size distribution that fits the scattered light images and SED of HH30 IRS (solid line). The dashed 
lines show ISM grain parameters (Kim, Martin, & Hendry 1994). The three panels show total opacity (upper), albedo (middle), and cosine 
asymmetry parameter (lower). 
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Fig. 2. — Evolution of SED with disk mass. Each panel shows model SED at ten viewing angles (evenly spaced in cosi), ranging from 
i = 87° (lowest curve) to i = 13° (upper curve). Also shown are the input stellar spectrum and optically thick flat and flared reprocessing 
isk models. The star is assumed to he at a. distance of SQ Opc. The short horizontal lines are 5cr SIRTF sensitivity limits for a 500 s exposure 



(http://sirtf.jpl.nasa.gov/SSC/B_Observing/SSC.B2.htm:). 
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Fig. 3. — Evolution of SED with disk mass. Each panel shows model SEDs for a fixed viewing angle for the range of disk masses, 
A^disk = 1O~ 8 M0 (lowest curve) to M^isk = 10 _1 Mq (upper curve). Input spectrum and SIRTF sensitivities are as in Fig. 3. 
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Fig. 5. — Color-color diagrams for our models. For each model we show twenty inclinations, evenly spaced in cosi. The numbers are the 
locations of different mass disks viewed pole on, 10 _1 Mq = 1, 10 _2 Mq = 2, 10~ 3 Mq = 3, etc. The very red colors are the highly inclined 
disks. As with Fig. 4, the bottom panel shows that long wavelength colors can distinguish a wide range of disk masses. 
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INFRARED SIGNATURES OF PROTOPLANETARY DISK EVOLUTION 

Kenneth Wood 1 ' 2 , C.J. Lada 2 , J.E. Bjorkman 3 , Scott J. Kenyon 2 , Barbara Whitney 4 , 

Michael J. Wolff 4 
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ABSTRACT 

We investigate the observational signatures of a straightforward evolutionary scenario for protoplan- 
etary disks in which the disk mass of small (<; 50/<m) particles decreases homologously with time, but 
the disk structure and stellar parameters do not change. Our goal is to identify optimal infrared spectral 
indicators of the existence of disks, their structure, and mass evolution that may be tested with the 
upcoming SIRTF mission. We present simulated spectral energy distributions and colors over a wide 
range of masses, 10~ 8 M Q < M^\sk < 10 _1 M©. Our Monte Carlo radiative equilibrium techniques enable 
us to explore the wide range of optical depths of these disks and incorporate multiple, anisotropic dust 
scattering. The SED is most sensitive to disk mass in the far-IR and longer wavelengths, which is already 
known from millimeter and radio observations. As the disk mass decreases, the excess emission of the 
disk over the stellar photosphere diminishes more rapidly at the longest rather than at short wavelengths. 
At near- infrared wavelengths, the disk remains optically thick to stellar radiation over a wide range of 
disk mass, resulting in a slower decline in the SED in this spectral regime. Therefore, near-IR excesses 
(K — L) provide a robust means of detecting disks in star clusters down to M^isk ~ 10~ 7 M Q , while the 
far-IR excess probes the disk mass, the caveat being that large inner disk holes can decrease the near-IR 
disk emission. 

Varying other disk parameters (outer radius, flaring, dust size distribution) alter the SED quantita- 
tively, but do not change our general conclusions on the evolution of SEDs and colors with the mass 
of small particles in the disk. Reducing the disk mass results in a clear progression in color-color 
diagrams with low mass disks displaying the bluest colors. We interpret color-color diagrams for Taurus- 
Auriga sources in the context of decreasing disk mass. Different viewing angles yield degeneracies in 
the color-mass relationship, but highly inclined disks are very faint and red and are readily identified in 
color-magnitude diagrams. 

Subject headings: radiative transfer — scattering — accretion, accretion disks — ISM: dust, extinction 
— stars: pre-main-sequence 



1. INTRODUCTION 

Protoplanetary disks are primarily detected through 
their signature infrared excess emission relative to a stel- 
lar photosphere (e.g., Mendoza 1968; Rydgren, Strom, & 
Strom 1976; Cohen & Kuhi 1979; Rucinski 1985). The IR 
excesses arise from the reprocessing of stellar radiation and 
the liberation of accretion luminosity within the disk (e.g., 
Lynden-Bell & Pringle 1974; Adams & Shu 1986; Adams, 
Lada, & Shu 1987; Kenyon & Hartmann 1987; Adams, 
Emerson, & Fuller 1990; Lada & Adams 1992; Hillenbrand 
et al. 1992). At other wavelengths, disk masses and ve- 
locity structures have been measured with radio and mm 
interferometers (Beckwith et al. 1990; Beckwith & Sar- 
gent 1993; Koerner, Sargent, & Beckwith 1993; Koerner 
& Sargent 1995; Dutery et al. 1996; Wilner & Lay 2000) 
and more recently, high resolution HST observations and 
ground based adaptive optics and speckle imaging tech- 
niques have imaged disks via scattered light (O'Dell, Wen, 
& Hu 1993; Burrows et al. 1996; Stapelfeldt et al. 1998; 
Padgett et al. 1999; Krist et al. 2000; Cotera et al. 2001; 
Grady et al. 2000; Roddier et al. 1996; Koresko 1998). 



The disks that are detected range in mass, structure, 
and formation mechanisms. Massive, optically thick, flar- 
ing circumstellar disks, M^\sk ^ 10~ 3 M Q , are a natural 
product of the star formation process (e.g., Shu, Adams, 
& Lizano 1987) and are detected around pre-main se- 
quence Classical T-Tauri stars. Low mass "debris disks" 
(Mdisk <, 1O~ 5 M ) are detected around older main se- 
quence stars (Backman & Paresce 1993). The dust that 
gives rise to the IR excesses and scattered light images of 
debris disks is believed to form from collisions of planetes- 
simals within the disks. 

The advent of large IR detector arrays now allows for 
simultaneous multiwavelength observations of large num- 
bers of disks in star clusters (e.g., Haisch, Lada, & Lada 
2000). The upcoming SIRTF mission will greatly increase 
the wavelength coverage for studying clusters, thereby en- 
abling detailed studies of disks over a wide range in age 
and mass. We therefore need to know what spectral indi- 
cators are best suited for detecting disks and determining 
their mass evolution. 

In this paper our working hypothesis is that through 
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time the disk mass decreases but the disk structure and 
stellar luminosity do not change. In our models the dom- 
inant opacity source for A <; 100/um is small (<; 50/um) 
grains, so another way of stating our model is that we as- 
sume the small particle mass traces the total mass of gas + 
dust. The small grain mass may decrease by the accretion 
of disk material or coagulation and growth of large grains 
and rocks while the gas is either depleted or accreted onto 
the star keeping the gas/dust ratio constant. The recent 
detection of i?2 in f3 Pictoris' disk and the derived gas/dust 
~ 100 is consistent with our assumption (Thi et al. 2001). 

We investigate the observational signatures of this sim- 
ple model for disk evolution and present radiative equilib- 
rium models for a range of disk masses and viewing angles. 
Model spectra are compared with semi-analytic approxi- 
mations for optically thick flat and flared disks and pre- 
dictions are made as to what spectral regions are most 
sensitive to the evolution of disk mass and at what wave- 
lengths disk detections are robust. We perform simula- 
tions for a fiducial disk structure and dust size distribution 
that fits the SED of the edge-on disk system HH30 IRS. 
(Wood et al. 2001). § 2 describes the ingredients of our 
models, § 3 presents our SED simulations and shows how 
SIRTF colors are sensitive to disk mass. § 4 presents mod- 
els that explore variations in parameters other than disk 
mass and discusses alternative disk evolution models. § 5 
compares our models with color-color diagrams of Taurus- 
Auriga sources and we summarize our findings in § 6. 

2. MODELS 
2.1. Radiative Equilibrium Calculation 

Theoretical temperature distributions and emergent 
spectra have been calculated for T Tauri disks using a 
variety of techniques. These include approximations for 
opaque flat disks (Adams & Shu 1986; Adams, Lada, & 
Shu 1987; Lada & Adams 1992) and extensions of this 
approach to flared disks (Kenyon & Hartmann 1987; Chi- 
ang & Goldreich 1997, 1999); two dimensional radiation 
transfer techniques (Efstathiou & Rowan- Robinson 1991); 
diffusion approximations (e.g., Sonnhalter, Preibisch, & 
Yorke 1995; Boss & Yorke 1996); disks which are modeled 
as "spherical sectors" (Men'shchikov & Henning 1997); 
vertical structure calculations performed by dividing the 
disk into plane parallel annuli (Calvet et al. 1992; Bell 
et al. 1997; Bell 1999; D'Alessio et al. 1998, 1999); 
and Monte Carlo techniques (Wolf, Henning, & Stecklum 
1999). Common approximations in many of the "tradi- 
tional" codes are to conduct the radiation transfer in only 
one direction in the disk and assume the dust grains scatter 
radiation isotropically. Also many techniques are limited 
to the study of optically thick and hence massive disks. 

Monte Carlo techniques are straightforward to adapt to 
any geometry, mass, and can accurately include polar- 
ization and multiple, anisotropic scattering. We use the 
Monte Carlo radiative equilibrium technique of Bjorkman 
& Wood (2001) which conserves energy exactly. For very 
optically thick disks we use our Monte Carlo technique for 
the upper layers of the disk (i.e., the disk "atmosphere") 
and a diffusion approximation for the densest midplane 
regions. For the most massive disks we simulate, typically 
less than 1% of the flux requires a diffusion treatment. 
More details of our adaptation of the Bjorkman & Wood 



(2001) technique for simulating T Tauri disks will be pre- 
sented in Bjorkman, Wood, & Whitney (2001, in prepa- 
ration). An advantage of using Monte Carlo techniques 
for studying a range of disk masses, is that we are not 
restricted to one-dimensional radiation transfer and can 
therefore simulate SEDs of low mass optically thin disks in 
which radial transport of photons is important. The out- 
put of our code is the disk temperature structure (due to 
heating by stellar photons and accretion luminosity) and 
the emergent SED and polarization spectrum at a range 
of viewing angles. A calculation of the hydrostatic disk 
structure (e.g., Chiang & Goldreich 1997; D'Alessio et al. 
1999) can be included in the Monte Carlo technique, but 
this would require an iterative scheme. At present we have 
not implemented such a scheme and instead perform the 
radiative equilibrium calculation for a fixed disk geometry. 

2.2. Disk Structure 

The determination of disk structure from fitting SEDs 
and scattered light images does not yield a single structure 
and dust size distribution that applies to all disks. Some 
systems are fit with passive flat disks (e.g., Adams et al. 
1987; Adams et al. 1990; Miyake & Nakagwa 1995), while 
others require flared disks heated by starlight and accre- 
tion luminosity (e.g., Kenyon & Hartmann 1987; Burrows 
et al. 1996; Stapelfeldt et al. 1998). The scattered light 
image of the edge-on disk of HH30 IRS (Burrows et al. 
1996) has for the first time allowed the vertical structure 
of a protoplanetary disk to be studied directly. We there- 
fore adopt the HH30 IRS disk as our fiducial model for our 
SED models. A fixed disk density structure which fits the 
scattered light images (Burrows et al. 1996; Cotera et al. 
2001) and SED (Stapelfeldt & Monetti 1999; Wood et al. 
2001) of HH30 IRS is 




where w is the radial coordinate in the disk midplane and 
the scale height increases with radius, h = ho (vd/R*)^ . 
For the HH30 IRS disk we adopt /3 = 1.25, a = 2.25, and 
h = 0.017.R*, giving /i(100AU) = 17AU. 

In our simulations the inner edge of the disk is trun- 
cated at the dust destruction radius, i?dust- Assuming 
T* = 4000K, R* = 2i? Q , and circumstellar dust subli- 
mates at 1600K, then i?dust ~ 8.R*. This is larger than 
the dust destruction radius for optically thin dust because 
the reprocessed emission from the disk provides additional 
heating over and above the direct stellar radiation, increas- 
ing the size of the dust destruction zone (a detailed dis- 
cussion of the shape of the dust destruction region will 
be presented in Bjorkman et al. 2001, in preparation). 
In currently popular magnetic accretion models the disk 
is truncated at a radius Rq, which may not be equal to 
^dust- If < -Rdust there will be a gas disk extending 
from _R to i?dust which may give rise to additional IR 
emission. Our models therefore assume that any material 
within i?dust is optically thin, which is a good approxima- 
tion for low mass disks. For high mass disks the gas may 
be optically thick producing larger near-IR excesses. 

2.3. Adopted Circumstellar Dust Properties 



Recent modeling of HST images (Cotera et al. 2001) 
and the SED of HH30 IRS (Wood et al. 2001) indi- 
cates that the circumstellar dust size distribution extends 
to larger grain radii than typical ISM grains. This is in 
agreement with many other observations indicating grain 
growth within protoplanetary disks (e.g., Beckwith et al. 
1990; Beckwith & Sargent 1991). This paper primarily in- 
vestigates the effects of disk mass on the SED and adopts 
circumstellar dust properties that reproduce the HH30 IRS 
SED. The dust model (chemical composition, mathemat- 
ical form for the size distribution, calculation of opacity 
and scattering parameters, etc) is described in Wood et 
al. (2001) and we only summarize the main features here. 
Specifically, we adopt, a size distribution 



n(a)da = Cia p exp (— [a/a c ] q ) da 



(2) 



with p = 3.5, q = 0.6, a c = 50/<m, a m ; n = 0.01/<m, and 
flmax = 1mm. The exponential scalelength, a c , yields dust 
particle sizes extending up to and in excess of 50//m. Fig- 
ure 1 shows the wavelength dependence of the opacity, 
scattering albedo, and Heyney-Greenstein phase function 
asymmetry parameter (Heyney & Greenstein 1941) for this 
size distribution. 

Recent Infrared Space Observatory spectra of Herbig 
Ae/Be stars (e.g., Meuss et al. 2000; Chiang et al. 2001; 
van den Ancker et al. 2000) are now allowing the circum- 
stellar dust chemistry to be studied. Because we adopt the 
dust properties of Fig. 1 for our simulations, we have not 
investigated the effects of different chemical compositions 
on the the resulting disk SEDs. How the circumstellar 
chemistry effects the SED is an interesting problem, but it 
is beyond the scope of this paper and we present models 
for different disk masses and the dust properties shown in 
Figure 1. In addition, our models do not include additional 
heating from transiently heated small grains (see the new 
radiation transfer code of Misselt et al. 2001) which is not 
important in the cooler Classical T Tauri stars considered 
in this paper. 

2.4. Energy Sources 

The energy input to the disk is from stellar photons and 
accretion luminosity liberated in the disk. As discussed 
in the previous section, we fix the disk structure for our 
radiation transfer simulations. Given the disk structure 
(Eq. 1), a disk theory determines the accretion rate for a 
given disk mass. Our parameterization of the disk den- 
sity and accretion follows that presented in the review by 
Bjorkman (1997), apart from the term ±^jR a /w. The 
accretion rate and viscosity parameter, M and ayisk, are 
related to the disk parameters by 



M = v / 18^3 a disk V c po hl/R* 



(3) 



where the critical velocity V c = \JGM±/R±. The flux 
due to viscous disk accretion, GM+M /2R+, is gener- 
ated throughout the disk midplane region according to 
(Shakura & Sunyaev 1973, Lynden-Bell & Pringle 1974) 



dE 
dAdt 



3GM+M 
4ttw 3 



(4) 



For low mass disks, the heating due to accretion luminos- 
ity is negligible and stellar irradiation dominates the disk 
heating. In general we choose aaisk = 0.01 (Hartmann et 
al. 1998), but for the most massive disk we simulate this 
results in a very large accretion luminosity, La > 0.8L*. 
For this case a d isk is adjusted so that La < 0.2L* in 
line with recent observational determinations of accretion 
luminosities in classical T Tauri stars (Hartmann et al. 
1998). 

3. MODEL RESULTS: SEDS, COLORS 

The following models use T* = 4000K, i?* = 2i? Q , 
-Rdisk = 100AU, and a distance of 500pc to the system. 

3.1. SED Evolution with Disk Mass 

Figure 2 shows the effect on the SED of changing the 
disk mass, but keeping the disk structure fixed to that of 
our fiducial model. The massive optically thick disks pro- 
duce SEDs that resemble that of the Kenyon & Hartmann 
(1987) flared disk model, aside from differences (addition 
of scattered light to pole-on views and silicate features) 
due to our inclusion of a finite albedo and non-gray opac- 
ity. The most massive disk has a very large near-IR excess 
due to the large accretion luminosity present in this model. 
The model SEDs display the characteristic features present 
in other simulations: large infrared excess emission, flat 
spectrum sources at intermediate inclinations, and double 
peaked spectra (optical and far-IR peaks) for very large 
inclinations. At large inclinations, the optical peak is due 
to scattered starlight as the dense disk totally obscures the 
star and disk emission at short wavelengths. 

Reducing the disk mass has the most dramatic effect 
at long wavelengths with the SED rapidly declining with 
decreasing mass. At short wavelengths the disk remains 
optically thick to stellar photons over a wide range of mass 
so the near-IR excess is not as sensitive to mass. Figure 4 
presents our models again, with the three panels show- 
ing SEDs for the range of disk masses at a given inclina- 
tion. These results show similar features to the models of 
Men'shchikov & Henning (1997, their Fig. 12) who pre- 
sented SEDs for a range of optical depths in a spherical 
geometry with evacuated bipolar cones. 

The SIRTF sensitivity limits show that at a distance 
of 500pc very low mass disks (M^isk J> 10~ 6 Mq) are de- 
tectable out to a wavelength of 70/<m. At 25/<m SIRTF 
will be sensitive to photospheric flux levels at 500pc, thus 
allowing for the detection of even lower mass disks. Our 
overall conclusion from these simulations is that near-IR 
excesses detect disks while far-IR excesses can be used to 
study their mass. 

3.2. Color Evolution with Disk Mass 

Figure 4 shows the variation of colors with disk mass 
for i < 60°. At long wavelengths we use the simu- 
lated flux at 70/um and 160/mi in forming the color and 
have not adopted any particular color system. The col- 
ors are fairly insensitive to inclination for i < 60°, but 
see §3.3 for color-color and color-magnitude diagrams that 
include highly inclined disks. For a passive disk with 
M disk = 1O -1 M , A(K - L) ?s 0.7, decreasing to a fairly 
constant A(K-L) ?s 0.4 for 10" 7 M Q < M disk < 10- 3 M Q . 
Larger K — L colors arise in massive disks where accretion 



luminosity is included, A(K—L) ps 1 for Mdi S k = 10 _1 M Q . 
In our models, less massive disks do not sustain large ac- 
cretion rates (see §2.3) and the K — L excess is due to 
reprocessing of starlight. For these disks, K — L is fairly in- 
sensitive to mass as the disk remains optically thick in the 
ncar-IR for masses M^\sk ^ 1O _7 M . Therefore ground 
based near-IR observations are capable of detecting very 
low mass disks. Again, we emphasize that this is the mass 
of small (<; 50/um) particles which dominate the near-IR 
opacity. 

At longer wavelengths the SED decreases with decreas- 
ing disk mass and this is reflected in the other color indexes 
remaining relatively flat with disk mass until the disk be- 
comes optically thin at the waveband under study. At the 
longest MIPS wavelength, 160/mi, the K — 160 color shows 
a clear progression from the most massive to least massive 
disk. 

3.3. Color-Color and Color- Magnitude Diagrams 

Figure 5 shows color-color diagrams for our model disks 
assuming intrinsic stellar colors from Kenyon & Hartmann 
(1995, Table A5). When we include all disk inclinations 
there are degeneracies in the color-mass parameter space, 
with large viewing angles generally leading to redder col- 
ors. This result differs from other investigations (e.g., 
Lada & Adams 1992; Kenyon, Yi, & Hartmann 1996; 
Meyer, Calvet, & Hillenbrand 1997) which assumed the 
disk emission was proportional to cos«, giving the bluest 
colors for edge-on viewing. More detailed radiation trans- 
fer modeling shows that a simple cos i scaling does not give 
the correct inclination dependence for flared disk models. 

The effect of changing the disk mass is not very strong 
in the (K - L)/(K - N) diagram for M disk £ 1(T 6 M , 
reflecting that at short wavelengths the disks remain op- 
tically thick over a wide range of mass. The disk masses 
separate out more clearly in the (K — N)/(K — 70) dia- 
gram. Therefore color-color diagrams of large numbers of 
sources that compare near and far infrared colors provide 
a means of determining the range of disk masses within a 
cluster. 

Figure 5 showed that inclination effects lead to a large 
spread in the location of different disk masses in the color- 
color diagrams. This leads to degeneracies in the lo- 
cation of different disk masses. However, Fig. 6 shows 
how color-magnitude diagrams can help in breaking the 
mass-inclination degeneracy. Placing all our disk models 
on color-magnitude diagrams shows that disks viewed at 
large inclinations occupy the lower righthand corners, i.e., 
highly inclined disks are faint and red. 

4. SED AND COLOR DEPENDENCE ON OTHER DISK 
PROPERTIES 

4.1. Flat Disks 

The previous sections presented SEDs for disks which 
have the flaring parameters that fit HH30 IRS images and 
SED. However, other investigations indicate a variety of 
disk structures in Classical T Tauri stars (e.g., Kenyon & 
Hartmann 1987; Miyake & Nakagawa 1995; Chiang et al. 
2001). To investigate the effects of flatter disks on SEDs 
and colors we repeated the simulations of Fig. 2 with a re- 
duced scaleheight h = 0.003i?* giving /i(100AU) = 3AU. 
The smaller scaleheight results in the disk intercepting less 



stellar radiation (e.g., Kenyon & Hartmann 1987) and con- 
sequently smaller excesses and bluer colors than our mod- 
els in § 3. For a given disk mass, the mm flux is unal- 
tered, but the mid to far-IR emission is sensitive to the 
disk structure (Kenyon & Hartmann 1987; Chiang et al. 
2001). The qualitative variation of the SEDs and colors 
with disk mass remains the same. 

4.2. Grain Growth 

In our models the dust size distribution does not change 
as the disk mass decreases. However, small dust grains 
may coagulate to form larger grains and rocks thereby al- 
tering the grain size distribution in the disk (e.g., Beck- 
with, Henning, & Nakagawa 2000). This will lead to 
changes in the wavelength dependence of the opacity. 
D'Alessio et al. (2001) investigated grain growth by keep- 
ing the disk mass constant and increasing the maximum 
grain size (effectively reducing the population of small 
grains) and changing the slope of the power law size dis- 
tribution. For a slope p = 3.5, their SED models show the 
largest effects occur for A ^ 20/<m, while forp = 2.5 (which 
puts more mass into larger particles) the SEDs are similar 
to ours in which the total disk mass decreases. SEDs aris- 
ing from either a low mass disk or a massive disk with large 
grains may be distinguished with detailed SED modeling. 

4.3. Disk Radius 

Decreasing the disk radius, but keeping the mass con- 
stant results in a disk that is denser and optically thicker 
than our models of § 3. A smaller disk effectively removes 
cool material at large radii that provides the bulk of the 
long wavelength emission, yielding smaller far-IR and mm 
fluxes (e.g., Beckwith et al. 1990). The other effect of 
squeezing the disk into a smaller volume is an increase in 
the height above the midplane at which stellar photons 
are absorbed. This gives a geometrically thicker dust disk 
that can intercept more stellar photons, raising the near 
and mid-IR emission. The overall effects of a very small 
disk are to increase the wavelength at which the disk be- 
comes optically thin yielding larger K — L colors and a 
slower decline in the K — 180 color with decreasing mass. 
The larger optical depth for small disks will lower the min- 
imum disk mass that may be detected from K — L colors. 
Large radius disks give the opposite of these effects: less 
near and mid-IR, and more far-IR and mm emission. 

4.4. Inner Disk Holes and Non-homologous Disk 
Evolution 

In order to investigate general trends, we have presented 
a simple model for disk evolution of protoplanetary disks 
in which the mass decreases homologously with time. Real 
disks will be more complicated than this, with the disk ra- 
dius growing, changes in the flaring parameters due to dust 
settling, and the possibility of accretion being terminated 
through the opening up of large inner disk holes and gaps, 
resulting in "inside-out" evolution. Current data does not 
provide much support for the existence of large holes in 
protoplanetary disks: if they do not show near-IR excesses, 
mid-IR and mm excesses are usually absent also (Stas- 
sun et al. 2001, Haisch, Lada & Lada 2001a), but there 
are some exceptions notably GM Aur (Koerner, Sargent, 
& Beckwith 1993). SI RTF should be able to distinguish 



between low mass protoplanetary disks, whose mass has 
evolved homologously (without creating large inner holes) , 
and disks which have evolved by creating large inner holes 
and thus directly test the inside-out disk clearing scenario. 

Debris disks often exhibit large inner holes or a ring-like 
structure (e.g., Koerner et al. 1998; Jayawardhana et al. 
1998; Schneider et al. 1999). However, debris disks are not 
likely to be protoplanetary disks that have evolved large in- 
ner holes. Debris disks are sufficiently old that the dust in 
them is not remnant protoplanetary material but instead is 
created and continually replenished by collisions between 
planetessimals formed previously in the disk. Moreover, 
protoplanetary disks which simply evolve from the inside 
out by clearing large inner holes (but otherwise maintain- 
ing their initial structure) , will still contain large reservoirs 
of material in their outer regions and thus should be con- 
siderably more massive than debris disks. 

5. COMPARISON TO OBSERVATIONS: DISK MASSES IN 
TAURUS-AURIGA 

Observations at mm and radio wavelengths, where the 
disk is optically thin, provide the best probe of disk mass 
(e.g., Beckwith et al. 1990). When such data is not avail- 
able our models show that far-IR, and in certain circum- 
stances even mid-IR data, can provide another means of 
probing disk mass (if the disks are not too massive and op- 
tically thick). For example, Kenyon & Hartmann (1995) 
presented a compilation of optical through far-IR obser- 
vations of Taurus- Auriga sources and noted a pronounced 
gap in the K — N distribution between the bluest Class II 
and the reddest Class III sources. Figure 7 shows the 
(K — L)/(K — N) color-color diagram for the Taurus- 
Auriga sources along with our face-on model colors for 
various disk masses. In the context of our models, the 
gap in the K — N distribution corresponds to disks with 
Md\sk <^ 1O~ 6 M0 and indicates that there are few circum- 
stellar disks with masses between 10~ 6 and 10~ 8 Mq in 
the Taurus-Auriga cloud (of course, this disk mass does 
not include rocks and planets). This, in turn, could be 
interpreted to indicate that once the disk mass falls be- 
low ~ 1O~ 6 M the timescale for clearing the remaining 
material in the inner (<; 5 — 10AU) disk is very rapid. 

6. SUMMARY 

We have investigated the observational signatures of an 
evolutionary model in which the disk mass decreases ho- 
mologously, but the disk structure and stellar parameters 
remain constant with time. Our main conclusions are: 

• Near-IR observations detect disks. Disks remain op- 
tically thick in the near-IR over a wide range of disk 
mass resulting in measurable K — L excesses for disks 
down to Mdisk ~ 1O _7 M . This corresponds to a 
dust mass of Mdust ~ 1O _9 M . 

• Mid-IR observations probe disk structure. Observa- 
tions in the 20/mi <; A <; 100/mi range are sensitive 



to the disk flaring parameters and are crucial for de- 
termining the degree of dust settling. 

• Far-IR most sensitive to disk mass. The far-IR emis- 
sion decreases rapidly with disk mass resulting in a 
strong correlation between K — 160 and disk mass. 
Our simulations indicate that at a distance of 500pc, 
SIRTF will be able to detect M disk £ 10" 6 M Q with 
a 500s exposure at 70/xm and even lower masses at 
25/im. 

• General trends are reproduced for a range of disk pa- 
rameters. Despite degeneracies in fitting SEDs of 
particular sources, we find that varying the disk pa- 
rameters within generally accepted limits introduces 
a spread in the color-mass relationship, but does not 
affect the general trends. 

The degeneracies inherent in fitting disk models to wave- 
length restricted datasets illustrates that it is necessary to 
have full SED coverage to constrain the disk structure and 
dust size distribution (e.g., Men'shchikov & Henning 1997; 
Chiang et al. 2001; D'Alessio et al. 2001). SIRTF, SMA, 
and ALMA will provide detailed SEDs for nearby sources 
enabling detailed modeling of the disk mass and struc- 
ture. For studies of stellar clusters, such data will not be 
available and mm observations, which are most sensitive 
to disk mass, cannot as yet achieve the required sensitiv- 
ity and resolution to study low mass disks and/or distant 
clusters. Therefore, we must appeal to statistics and de- 
termine average properties of sources in different clusters 
for studying disk evolution. SIRTF will provide broad 
band colors for hundreds of sources within a given cluster, 
from which we can construct median colors and the spread 
around the median for each cluster (e.g., Haisch, Lada, & 
Lada 2001b). From our models it is apparent that NIR ex- 
cesses detect disks, and therefore the absence of a K — L 
excess within a cluster of a given age can determine the 
maximum disk lifetime. By comparing the median colors 
from different clusters with a spread in age (as a whole 
and as functions of position within the cluster), we can 
search for trends as a function of cluster age. Therefore, 
the far-IR probed by SIRTF may allow us to determine 
a timescale for disk evolution — i.e., in our homologous 
evolutionary model we would be able to convert the color- 
mass correlation to one which tracks disk mass in clusters 
as a function of age, analogous to the recent ISO analysis 
of debris disks by Spangler et al. (2001). 
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Fig. 1. — Dust parameters for a grain size distribution that fits the scattered light images and SED of HH30 IRS (solid line). The dashed 
lines show ISM grain parameters (Kim, Martin, & Hendry 1994). The three panels show total opacity (upper), albedo (middle), and cosine 
asymmetry parameter (lower). 
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Fig. 2. — Evolution of SED with disk mass. Each panel shows model SED at ten viewing angles (evenly spaced in cosi), ranging from 
i = 87° (lowest curve) to i = 13° (upper curve). Also shown are the input stellar spectrum and optically thick flat and flared reprocessing 
disk models. The star is assumed to be at a distance of 500pc. The short horizontal lines are 5a SIRTF sensitivity limits for a 500 s exposure 
(http://sirtf.jpl.nasa.gov/SSC/B_Observing/SSC.B2.html). 
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Fig. 4. — Colors as a function of disk mass for inclinations i < 60°. Triangles are passive disks, squares include accretion luminosity. This 
bottom panels show that long wavelengths are more sensitive to disk mass. 





Fig. 5. — Color-color diagrams for our models. For each model we show twenty inclinations, evenly spaced in cosi. The numbers are the 
locations of different mass disks viewed pole on, 10 -1 Mq = 1, 10~ 2 Mq = 2, 10~ 3 Mq = 3, etc. The very red colors are the highly inclined 
disks. As with Fig. 4, the bottom panel shows that long wavelength colors can distinguish a wide range of disk masses. 
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Fig. 6. — Color-magnitude diagrams for our models assuming a distance of 500pc. The edge-on disk sources are very faint and red and 
occupy the lower right corners of the color-magnitude diagrams. 



3.0 



2.5 



2.0 



+ ++ 



X 



1 .5 



+ * ; 



1 .0 



0.5 



++ + ++ + + + 



+ + 



+ + + + * 2+ 

+ * + + + if + 

■m- + + 



0.0 



o* o 



K - N 



Fig. 7. — Color-color diagrams for Taurus- Auriga sources Class III (diamonds), Class II (crosses), Class I (stars). The data are taken from 
the Kenyon & Hartmann (1995) compilation. The numbers are our models for different mass disks viewed pole-on as in Fig. 5. The gap in 
the K — N distribution between the Class II and Class III sources is filled in by our models having M^k < 10~ 6 Mq, which may imply that 

the timescale to clear M<ji s k < 10 -6 Mq is very rapid. 



